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The anomalous low energy behaviour observed in metals with strong electron correlation, such 
as in the heavy lermion materials, is believed to arise from the scattering of the itinerant electrons 
with low energy spin fluctuations. In systems with magnetic impurities this scattering leads to the 
Kondo effect and a low energy renormalized energy scale, the Kondo temperature Tk. It has been 
generally assumed that these low energy scales can only be accessed by a non-perturbative approach 
due to the strength of the local inter-electron interactions. Here we show that it is possible to 
circumvent this difficulty by first suppressing the spin fluctuations with a large magnetic field. As 
a flrst step field-dependent renormalized parameters are calculated using standard perturbation 
theory. A renormalized perturbation theory is then used to calculate the renormalized parameters 
for a reduced magnetic field strength. The process can be repeated and the flow of the renormalized 
parameters continued to zero magnetic field. We illustrate the viability of this approach for the 
single impurity Anderson model. The results for the renormalized parameters, which fiow as a 
function of magnetic field, can be checked with those from numerical renormalization group and 
Bethe ansatz calculations. 



I. INTRODUCTION 

Materials which have been classified as having strong 
electron correlation are those which have itinerant elec- 
trons in narrow bands arising from atomic d and f states. 
The effective on-site inter-electron interaction for elec- 
trons in these shells is of the order, or greater than, the 
band widths so that the direct application of perturba- 
tion theory to models of these systems, such as the Hub- 
bard model or periodic Anderson model, is not valid in 
most parameter regimes — particularly in the regimes 
where anomalous behaviour is expected. The same ap- 
plies for impurity models, such as the Anderson impu- 
rity model, where the anomalous low energy scale, the 
Kondo temperature Tk, occurs only when the impurity 
electron is almost localised and experiences a strong on- 
site interaction. For this reason the search has been for 
non-perturbative techniques that can handle the strong 
on-site interactions. Such non-perturbative techniques as 
the numerical renormalization group (NRG)^'^ and Bethe 
ansatz (BA)-^'^ have had considerable success in tackling 
most impurity problems, and some significant progress 
has been made in understanding the behaviour of lat- 
tice models, such as the Mott-Hubbard metal-insulator 
transition, using dynamical mean field theory (DMFT)^ 
in conjunction with a non-perturbative impurity solver. 
Many other aspects of strong correlation behaviour, such 
as quantum criticality and the competition between vari- 
ous low temperature broken symmetry states, both mag- 
netic and superconducting, are still very much open prob- 
lems where there is a need for new ideas and techniques. 
A feature of strongly correlated systems is that there are 
low lying electronic collective excitations, which scatter 
strongly the itinerant electrons. In systems with a strong 
on-site interaction U{> 0), these are usually spin excita- 
tions, arising from the effective ferromagnetic or antifer- 
romagnetic exchange interactions, which appear explic- 
itly in models such as the t-J models, which can be de- 



rived from a Hubbard model when U is large, where the 
antiferromagnetic exchange term J ^ t'^/U and t is the 
hopping matrix element. 

Here we put to the test a new approach for accessing 
the low energy scales by applying it to an impurity An- 
derson model in the strong correlation regime. We have 
accurate results for this model from the Bethe ansatz 
(BA) and the numerical renormalisation group (NRG) 
which we can use to check our results. The motivation 
is to develop a technique which we can apply to a wider 
class of models. The Bethe ansatz, though powerful, is 
restricted to integrable models, which cover a wide class 
of impurity models, and some one dimensional models 
such as the Hubbard model, but there is no generali- 
sation to models of higher dimension. The NRG can 
only be applied to impurity models which have a low 
degeneracy of impurity states so that the matrix sizes 
do not become too large for practical iterative diagonal- 
isation. In conjunction with the DMFT the NRG can 
also be applied to infinite dimensional models which can 
be mapped into effective impurity models. The density 
matrix renormalisation group (DMRG)^ is an alternative 
form which has been applied successfully to one dimen- 
sional systems. There are also functional forms of the 
renormalisation group (fRG)^ and non-perturbative RG 
approaches which are at present being developed and ap- 
plied with some success, but the problem of dealing with 
the strong correlation of models for systems in two and 
three dimensions is still an open challenge. 

The approach we consider here is based on a renor- 
malised perturbation theory (RPT). The RPT has been 
developed and tested in detail for the Anderson impu- 
rity model-i^. We begin by giving a brief description of 
this application and then consider the problems in ap- 
plying the method more generally, and how they might 
be overcome. 
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The Hamiltonian for the Anderson model^'' is 

Ham = ^ ed,adl,d^ + U nd,•t■f^d,i (1) 

cr 



where ed.o- — ^d — <^g^^BH /2 is the energy of the locahsed 
level at an impurity site in a magnetic field H , U the 
interaction at this local site, and Vk^a the hybridization 
matrix element to a band of conduction electrons of spin 
a with energy ek,a = Cfe — o^gc^J■BH/2, where gc is the 
g-factor for the conduction electrons. When U = the 
local level broadens into a resonance, corresponding to 
a localised quasi-bound state, whose width depends on 
the quantity A^{uj) = TrJ2k \Vk,cr\'^5{uj - ek,a)- For the 
impurity model, where we are interested in universal fea- 
tures, it is usual to take a wide conduction band with a 
flat density of states so that Act(w) becomes independent 
of a;, and can be taken as a constant Ao-. In this wide 
band limit Acr(w) will be independent of the magnetic 
field on the conduction electrons, so we can effectively 
put Qc — Q. When this is the case Ao- is usually taken 
to be a constant A independent of a. We also introduce 
the notation h — gpi^H /2. 

In the renormalized perturbation theory approach^ 
we cast the corresponding Lagrangian for this model 
J^AM{^d,(T, A, U) into the form, 

^AM(ed.(T, A, U) = £AM(ed, cr, Act, [/) -f £ct (Al,cr , A2, cr, A3), 

(2) 

where the renormalized parameters, Id^a and Ao-, are 
defined in terms of the self-energy ll„{uj,h) of the one- 
electron Green's function for the impurity state, 
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(3) 



- eda + iA - Scr(w, h) ' 

and are given by 

ed,a ^ zaied^a + ^aiO,h)), ^ z^A, (4) 

where is given by z^r — 1/(1 — h)). The renor- 

malized or quasiparticle interaction U{h), is defined in 

terms of the local total 4- vertex T^^^{uji,uj2t'^3,i^a] h) at 
zero frequency, 



C/(/i) -Zf^;r|J(0,0,0,0;/i). 



(5) 



It will be convenient to rewrite the spin dependent quasi- 
particle energies in the form, ed.cr — id{h)~ahf}{h), where 

^d{h) = ed,<T, f}{h) = cr?d,cr, (6) 



where ed{h) and f}{h) are both even functions of the mag- 
netic field h. 



In terms of the renormalised parameters, the Green's 
function takes the form, 



uj — ed{h) + ah + iAcr{h) — Scr('^, h) 



(7) 



where h = f)(h)h, and T,„{uj^h) is the renormalised self- 
energy given by 

t,{Lo, h) = z,(S,(c^, h) - I],(0, h) - c^E;(0, h)). (8) 

The propagator in the RPT is the free quasiparticle 
Green's function. 



G.(a;,/i) 
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ijj — ed{h) + (rh + iA^{h) 



(9) 



The expansion is carried out in powers of the interac- 
tion U{h) for the complete Lagrangian defined in equa- 
tion The counter term part of the Lagrangian 
>Cct(Ai^cr, A2,CT, A3), given by 

>Cct(Al,.,A2,a,A3) = (10) 
'^cr(T)(A2,CT'^ + >^l,a)dcr{T) + A3nd,-f (t)^^,^ (t) , 



E 



where da{T)dc,{T) = nd.a-{T), and d„[T) and d„{T) are the 
Grassmann fields corresponding to the impurity creation 
and annihilation operators, which are integrated over in 
the calculation of the partition function. This part of 
the Lagrangian essentially takes care of any overcount- 
ing. The parameters, ed,cr, A^. and U , have been taken 
to be the fully renormalized ones, and the counter term 
parameters, Ai^a-, A2,ct and A3, are required to cancel any 
further renormalisation. This leads to the renormalisa- 
tion or over-counting conditions, 



Sc.(0,/i) -0, 



and 



rt,i(0,0,0,0;/i) = f/(/i). 



(11) 



(12) 



The counter terms are completely determined by these 
conditions. 

Exact results for the Fermi liquid regime, which 
were first derived in a phenomenological approach by 
Nozieres^^, can be derived from the RPT, working only to 
second order in U . We give some of these results for the 
particle-hole symmetric model. In this case, id{h) = 
and Ag. is independent of cr, so we drop the a index for 
this quantity. The free quasiparticle density of states 
Pa{^^^ h), given by the spectral density of the free quasi- 
particle Green's function in equation ([9]), takes the form. 



p,{uj,h)=p^^\uj,h,A{h)), 



(13) 
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where p''a^\uj,h, A) is the local density of states for the 
non-interacting system, 



A 



(14) 



As p^\lu, h, A) becomes independent of for a; = 0, we 
can drop the a index when w = 0. 

The impurity contribution to the coefficient of the elec- 
tronic specific heat 7(/i) is directly proportional to the 
quasiparticle density of states, 



(15) 



for u! 0, where Us{h) is the effective quasiparticle in- 
teraction in the longitudinal channel given by Us{h) = 
U{h)/{l + U{h)p{0,h)). 

The corresponding expression for the transverse sus- 
ceptibility ('^j ^) which we denote by Xt(<^) ^) is 



nt(w,/i,A(/i)) 



\-Ut{h)I].t{i^,h,A{h)) 



(23) 



where 



1 , fiA-h\ 
-In — oj = -2/i, 



7r(iA-/i) 27rA \iA + h 



where fee is the Boltzmann constant. The induced mag- 
netisation M{h) is given by M{h) = g^B'm{h), where 



m{h) = i((ndt) - M) = moik A{h)), (16) 

where {n^a-} is the expectation value of the occupation 
number of the impurity site, rido-, and mo(/i. A) is the 
magnetisation for the non-interacting model given by 



ma{h,A) = - tan ^ ( 

TT V A 



(17) 



The longitudinal susceptibility XiW (in units of (ff/ie)^) 
is given by 

Xiih) = p(0, h){l + U{h)p{Q, h))/2, (18) 
and the local charge susceptibility Xc{h) is given by 

Xc{h) = 2p{Q,h){l-U{h)p{Q,h)). (19) 
The perpendicular susceptibility Xi(^) is given by 

mih) 



2h ' 



(20) 



in the limit of zero transverse field^^. 

Asymptotically exact results can also be derived for 
low energy dynamic spin and charge susceptibilities by 
taking into account repeated quasiparticle scattering^^. 
The longitudinal dynamic spin susceptibility is given by 



yAuj, h) — — 

2{l~Us{h)niiLu,h,Aih))) 

for the symmetric model, where 



(21) 



Ili{uj,h,A) 



-A 



■7TUj(uj + 2iA) 



In 



nih^ + A^) 



uj + iA — h 
iA - h 



for w = 0, 



In 



uj + iA + h 
lA + h 



(22) 
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-In 



uj + iA + h 



TT \ bj + 2h + 2iA V lA + h 



1 ^^( u; + ^A + h ^^^ u.^-2h.{24} 



L0 + 2h 



iA - h 



The effective quasiparticle interaction in this channel 
Ut{h) is determined by the condition that the dynamic 
transverse susceptibility at w = it is equal to twice the 
static perpendicular susceptibility given in equation pop , 
Xt{0,h) = 2x±ih). This gives 



Xt{0,h) 



mo (/i, A) 
h 



mo{h,A) 



Ut{h)mo{h,A) 



(25) 



on using the result llt{0,h, A) = mo(/i, A)//i for cj = 0. 
We then find Ut{h) must be such that 



h = h + Ut{h)m{h). 



(26) 



With Ut{h) defined as in equation ([25l) . the exact expres- 
sion for the magnetisation takes the form of a mean field 
equation, 



1 



m{h) — —tan 

TT 



h + Ut{h)m{h) 



(27) 



The only equations that we have quoted that are not 
exact are the equations for the dynamic susceptibilities 
Xii^J^h) and Xt(w,/i) for finite frequency w 7^ 0. They 
have nevertheless been shown to provide a very good ap- 
proximation to the NRG results for these quantities over 
the whole low frequency range and satisfy the Korringa- 
Shiba relation exactly (for more details see referenceii^) . 

To evaluate these formulae we need to know the renor- 
malised parameters. In the Kondo regime, in the ab- 
sence of a magnetic field, these can be reduced to a single 
parameter—, the Kondo temperature Tk (xi(0) — ^/'^Tk) 
such that 



U{0) = 7rA(0) = 4Tk, 



(28) 



but this parameter is still undetermined. The most ac- 
curate way of calculating the parameters in terms of the 



'bare' parameters of the model, e^, A and C/, is an indi- 
rect one, from an analysis of the low energy fixed point 
in an NRG calculation. The procedure for doing this is 
described elsewherei^ . This approach works very well 
but restricts the method to models where we can apply 
the NRG, which means one where we already have a so- 
lution. It is, nevertheless, a useful adjunct to the NRG, 
enabling one to calculate many quantities more easily and 
more accurately. However, we want to develop a way of 
calculating the renormalised parameters which is inde- 
pendent of the NRG and can be applied more widely. 
The obvious method would be to calculate them directly 
from their definitions in equations (|4]) and (O, but this 
would again seem to require having a solution, at least 
for the low energy behaviour of the self-energy. A direct 
perturbation approach to calculate the self-energy would 
not seem to offer the possibility of accessing the strong 
correlation regime, as we know for the symmetric model, 
low order perturbation theory is unreliable for JJ/tt A > 1 , 
and no-one has so far succeeded in summing a subclass 
of terms to give the correct low energy behaviour. 

We do know, however, that for this model the large 
renormalisation effects arise for large positive U from 
the scattering of local low energy spin fluctuations. If 
we apply a strong local magnetic field these spin fluctu- 
ations are suppressed, so the renormalisation effects are 
much weaker, and perturbation theory can be applied in 
this case, even in the regime U/wA ^ 1—. This im- 
plies that in the limit of a very large magnetic field, we 
can evaluate the renormalised parameters directly from 
equations Q and (O, using perturbation results for the 
self-energy. This assumption can be checked, as we have 
results for the Anderson model of the renormalised pa- 
rameters for any value of the magnetic field, deduced us- 
ing the NRGi^"— . An example is shown in figure[l]for the 
symmetric Anderson model, with a value of U/nA = 3, 
ttA — 0.1, such that in very low magnetic fields we are 
in the strong correlation regime with a Kondo tempera- 
ture Tk = 0.002. The degree of renormalisation in zero 
field can be estimated from 1/z — A/A, corresponding 
to a mass enhancement factor, which in this case gives 
1/z ^ 12.7. We can see from figure [T] that z increases 

monotonically with increase in the value of the magnetic 
field. It does, however, need extremely large field val- 
ues, such that h ^ U, before the renormalisation effects 
are completely suppressed. The other two parameters, 
77(/i) and U{h), do not simply increase monotonically as 
the magnetic field is increased. Initially they increase 
rather slowly, and then rise much more rapidly to a peak 
value, and then in the extreme large field limit approach 
their bare values. The initial increase is not surprising 
because the spin fluctuations are being suppressed and 
the quasiparticles are becoming less renormalised, the 
first stage of their undressing. We see that the peak 
value of U{h) is greater than the bare value U. This 
is because at this point we are approaching the regime 




f) h -Jr--. J = ri----|= = [ T I , \ , I , \ , \ , \ , 
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ln(h/Tjj) 

FIG. 1: Plots of the renormalised parameters, A{h)/A, 
U{h)/U and fi{h), for the symmetric Anderson model, cal- 
culated using the NRG for the case ttA = 0.1, U/nA = 3 as 
a function of the logarithm of the magnetic field h/Tn, where 
Tk = 7rA(0)/4 = 0.002. 



where mean field theory and RPA are applicable. In this 
regime, RPA corresponds to substituting Us{h) — > U and 
Ut{h) — ^ J7, with A(/i) — >• A into equations ([2T|l and 
(j23p . giving h h — h + Um(h) and from equation 
(fT8|. U{h) U/{1 - C/p(°)(0,/i, A)). This can explain 
why U{h) can be enhanced over the bare value U . As 
p^°^(0,ft-, A) — > as ft, — 7> oo, the enhancement disap- 
pears for very large field values. The other term con- 
tributing to U{h) from the z-factor in ([5]) has the oppo- 
site effect. In the mean field regime for very large fields, 
where z{h) = 1, it has no effect but plays a dominant role 
in reducing U{h) in the weak field regime where z{h) is 
small. 

The peak in fi{h) can be explained in a similar way. For 
large fields we have the mean field result for fi{h) = 1 + 
Um{h)/h, so that fj{h) — >■ 1 as /i — oo because m{h) — 
1/2. It follows from equations and (fT5|) that in the 
limit /i — > the value of 77(0) is equal to the Wilson x/l 
ratio, which in the strongly correlated (Kondo) regime 
takes the value 2. The fact the z{h) increases as for 
small h explains the initial increase of fj{h) with h. That 
a peak should occur at an intermediate field value follows 
as the likely behaviour from the extrapolated trends at 
low and large field values. 

We have conjectured that we should be able to explain 
the NRG results for the renormalised parameters in the 
large magnetic field limit using the leading perturbational 
corrections arising from mean field theory and RPA. The 
question then arises: Can we find a way to continue the 
process to lower magnetic fields, as the renormalised pa- 
rameters in large fields are continuously connected to the 
strongly renormalised values in the weak field regime? 
We know that if we start at the other limit, with the 
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zero field renormalised parameters known, we can use the 
RPT to calculate the renormalised self-energy in a weak 
field h^. From this result, we could then calculate the 
change in the renormalised parameters due to the intro- 
duction of the weak magnetic field. We could generalise 
this idea by considering that we know the parameters 
for an arbitrary magnetic field value h, and use these 
to calculate the self-energy for a system with a slightly 
smaller or larger magnetic field h ± Sh, and hence de- 
duce the renormalised parameters for these neighbouring 
magnetic field values. This would give a set of scaling 
equations for the renormalised parameters as a function 
of the magnetic field. If we can obtain good starting val- 
ues from perturbation theory in the large field limit, we 
should be able to reduce iteratively the magnetic field, 
and scale into the strongly correlated weak field regime. 
The aim of this paper is to test the feasibility of this 
scheme, and we use the results from the NRG and Bethe 
ansatz to test any approximation. We begin first of all 
looking at the behaviour in the limit of large magnetic 
field. 



II. LARGE MAGNETIC FIELD LIMIT 



A diagrammatic representation of this equation is given 
in figure [5J The skeleton diagrammatic formulation of 
the perturbation theory is particularly useful when self- 
consistent approximations are used. The simplest self- 
consistent approach is the mean field theory where only 
the tadpole diagram, corresponding to the first term in 
equation is taken into account. For the particle-hole 
symmetric model this gives So-(a;) — -~aUm{h), where 
m{h) is determined self-consistently from the equation. 



m{h) 



1 



■ tan 



(30) 



where h = h + Um{h). Note that we have not included 
the constant term U/2 in the self-energy because for the 
symmetric model it can always be absorbed into the en- 
ergy level Ed {— —U/2) to give — 0. From equation 
this implies that in the mean field regime Ut{h) — ?> U. 
For this approximation we deduce the renormalised pa- 
rameters, fj{h) and A(/i). As the self-energy is indepen- 
dent of uj, it follows that z{h) — 1, so 



,(/.). 1 + -tan- - 



Hh) = A, 



(31) 



A formally exact expression can be written down for 
the self-energy Sct(i^) in terms of skeleton diagrams, 
where the propagators correspond to the full many-body 
Green's function Gd.ai'^), as given in equation ((S]), and 



the full vertex r[,.'*|'_o.(cji, cj2, 1^3, 1^4), 



duj' duj" 
27ri 27ri 



(29) 




CO CT 



CT CB+CO' CT 



(4) 

J~'[o)+co',ci)"-(ii',ci)",co] 



CO CT 



FIG. 2: The skeleton diagrams for the self-energy T^criuj) in 



where m{h) is determined self-consistently from (pO|) . We 
can deduce the static longitudinal susceptibility by dif- 
ferentiating m{h) with respect to ft,, and from equation 
(UH]) deduce U . The result is Us{h) — U, or equivalently, 



U{h) = 



u 



l-C/p(o)(0,ft,A)' 



(32) 



We anticipated this result in the previous section, based 
on the RPA approximation for the longitudinal dynami- 
cal susceptibility. 

We now have a 1-1 correspondence between the RPT 
results and the mean field/ RPA equations in the trans- 
verse spin scattering channel, and we conjecture that 
asymptotically for very large magnetic field values, 

h^h = h + Um{h), A{h)^A, Ut{h)^U. (33) 

In figure [3] we compare the mean field results for z{h) = 
A(/i)/A, Ut{h)/U and fi{h), with the corresponding re- 
sults calculated using the NRG as a function of ln(/i/rK), 
for the model for ttA = 0.1 and U /ttA = 3. As expected 
the mean field results are in agreement with the NRG re- 
sults for asymptotically large magnetic fields. The mean 
field results give a good approximation for these param- 
eters for values of ft > QOTk (ln(ft/rK) > 4.5). Though 
the mean field results are only valid for very large mag- 
netic field values, we can build upon them by expanding 
about the mean field solution. The propagators in the 
perturbation expansion now take into account the mean 
field self-energy so that the tadpole diagrams, or mean 
field insertions, no longer appear explicitly. 



terms of the 4-vertex F 



(4) 
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FIG. 3: A comparison of the renormalised parameters, f]{}i) 
z(h) — A{h)/A and Ut{h)/U, from mean field tlieory with 
those calculated using the NRG, for the symmetric Ander- 
son model, with ttA = 0.1, U/nA = 3 as a function of the 
logarithm of the magnetic field h/T-K, where Tk = 7rA(0)/4. 



For U > ttA in the absence of a magnetic field, the 
mean field solution predicts a state with a local mag- 
netic moment, such that it costs no energy to flip the lo- 
cal moment. This is reflected in the dynamic transverse 
spin susceptibility, calculated with mean field propaga- 
tors, which develops a singularity at w = 0. Hence, the 
most important corrections to the self-energy are likely to 
arise from these spin flip scattering processes. In evaluat- 
ing the self-energy using equation (|29p we need to include 
the dominant spin flip scattering terms contributing to 
the 4- vertex t''^^_^{uj + uj' , uj" — uj' , uj" , uj) , which are illus- 
trated in figure |31 These diagrams are the same as those 
which are taken into account in the RPA expression for 
the transverse dynamic susceptibility Xt{^)- Summing 
this class of diagrams gives as an approximation for the 



4- vertex, r!j,_|(w + cj', cj" - cj', w", w). 



U 



1 - Unt{uj',h,A) 



(34) 



Substituting this result into equation (f29| gives the re- 
sult, 



duj' 
2^' 



(35) 



where the Green's function G'™^(aj,/i) is the mean field 
propagator and the transverse dynamic susceptibility 
X((cj, h) is calculated in the RPA. There is a similar ex- 
pression for the spin down self-energy, but with particle- 
hole symmetry they can be reduced to a single equation 
as G^(w, h) = -G|(-w, h). 

As the self-energy acquires a dependence on the fre- 
quency uj in this approximation, in the calculation of 
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FIG. 4: The RPA diagrams that contribute to the 4-vertex 
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the renormalised parameters, there will be a change of 
the quasiparticle weight factor z(h) from the mean field 
value 1. We can see from the results for the renormalised 
parameters calculated using this result, shown in figure 
[51 that this leads to a significant improvement in com- 
paring the results with those calculated using the NRG. 
The renormalised parameters are plotted as a function 
of hi{h/T]^) for the same parameter set as in figure |31 
ttA = 0.1, U/nA = 3. The values for both z{h) and 
Ut{h) are now seen to decrease with a decrease in ft, in a 
rather similar way to the NRG results. The value of fi{h) 
also develops a peak as the magnetic field value is reduced 
in a similar way to the NRG results, but the calculations 
for this quantity break down for ln(/i/TK) < —0.5. The 
results now constitute a good approximation to the NRG 
results down to a magnetic value h ^ 27T]^ corresponding 
to ln(/i/Tk) ~ 3.4. 



-- z(h)RPA 

ri(h) RPA 

- .- D|(h)/nARPA 
-- z(h)NRG 

TJCh) NRG 

.,_ U(h)/jtANRG 




FIG. 5: A comparison of the renormalised parameters, fi{h) 
z(h) — A{h)/A and Ut{h)/U, using the RPA propagator in 
calculating the self-energy with those calculated using the 
NRG, for the symmetric Anderson model, with irA = 0.1, 
U/ttA = 3 as a function of the logarithm of the magnetic 
field h/TK, where Tk = 7rA(0)/4. 
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We have derived, using perturbation theory, asymptot- 
ically exact expressions for the renormalised parameters 
for very large magnetic field values. However, the mag- 
netic field regime where these results are valid is com- 
pletely outside the range that can be realised experimen- 
tally, except possibly for some quantum dot systems. We 
need a way of extending the calculations to much lower 
magnetic field values. It should be possible, in princi- 
ple at least, to improve on these results by systemati- 
cally taking corrections to the RPA into account. How- 
ever, this gets more and more difficult, involving multiple 
frequency integrations and more complex self-consistent 
equations to solve. Our aim here is a limited one, the 
calculation of the renormalised parameters which depend 
only on the form of the self-energy in the low frequency 
regime, and as a result we can adopt a different strategy 
for carrying out higher order calculations. We return to 
the idea discussed at the end of the Introduction, that we 
might be able to use the RPT, with renormalised param- 
eters for a magnetic field h, to calculate the renormalised 
self-energy at a reduced fields value h — dh, and hence de- 
rive the renormalised parameters for the lower magnetic 
field value h — Sh. 



III. CALCULATIONS USING THE RPT 

The RPT approach differs from standard perturbation 
theory as it includes counter terms. We should con- 
sider how to handle these terms, and test various ap- 
proximations for calculating the renormalised self-energy 
Yi^{u},h) in a given magnetic field h, using the renor- 
malised parameters, fi{h), A(/i) and U(h), for the same 
magnetic field value before considering how to use 
them for a system with a reduced magnetic field. 

When RPT calculations are carried out to a specific or- 
der n in powers of {/, we can handle the counter terms by 
also expanding them formally in powers of J7, and taking 
into account all diagrams generated to order n, including 
those involving the counter terms. Counter term contri- 
butions of order greater then n will not contribute and 
the terms up to order n can be determined, order by or- 
der, by requiring them to satisfy the renormalisation con- 
ditions given in equations pT|) and ((T2|) (see reference^ for 
a specific example). We need a more general procedure, 
however, for calculations based on summing subclasses of 
diagrams taken to infinite order. 

There are five counter terms to deal with, four of them, 
Ai.o- and \2.a involve one-body terms, and the fifth is an 
interaction term A3. The contributions arising from the 
diagrams for the counter terms, Ai.o- and A2,ct, alone, can 
be taken into account fully as they do not involve the 
interaction term. The interaction counter term A3 can 
be added to the interaction U , and the expansion carried 
out in powers of the total interaction term J7 -I- A3 . The 
value of A3 is then required to satisfy the renormalisa- 
tion condition in equation (IT^ . The local Green's func- 
tion subject to a magnetic field ft,, using the renormalised 



parameters for the field /i, takes the form. 



(36) 



z(h) 



oj + ah + iA{h) + Ai,„ + cjAs,^ - S^(w, h) 



where E^(a;,/i) is the self-energy calculated using the 
quasiparticle propagator, Ga{i^,h) = Gr,{uJ,h) / z(h), 
which includes the two counter terms, Ai ,j and A2,(t. This 
means that the self-energy Sf (1^, h) will be a function of 
Ai^CT and as well as A3, so that the renormalisation 
conditions in equations take the form. 



S^(0, /l, Ai^o-', A2,ct', A3) — Ai^cr, 



and 



5f;f(i:j,/i, Ai,o^', A2,<T',A3) 



D = 



(37) 



(38) 



Given a result for Sf(0, ft), calculated from a particu- 
lar subset of diagrams, the conditions in equations ([37]) 
and (pS)) generate self-consistent equations which must 
be solved to determine the counter terms, Ai.o- and A2,o-. 

We can now perform calculations equivalent to mean 
field and the RPA but in the RPT framework. At the 
first mean field stage we only take the tadpole diagrams 
into account giving S^(a;,ft) = —(T{U + X'i)m{h). In 
this approximation from equations p7[) and psp we get 
Ai,CT = cr(C/ -I- A3)m(ft) and A2,cr = 0, so 



GTi^.h) 



1 



crft + iA(ft) 



(39) 



The next stage is to calculate the transverse susceptibility 
Xt(w, ft), using propagators which include the mean field 
insertions. The Ai.o- in the propagator given in equation 
(|57p cancels the mean field term, the propagator now 
becomes the free propagator given in equation The 
result for Xt{'^^ h) is 



Xt(w,ft) 



nt(w,ft,A(ft)) 



l-(C/ + A3)nt(c^,ft,A(ft)) 



(40) 



The counter term A3 has yet to be determined. We have 
the exact result, xt(0, ft) = m(ft)/ft, which from (j40|) im- 
plies ft + (C/ -I- A3)m(ft), so we can identify [/ -I- A3 as 
Ut{h) in equation (j26p and the expression for Xt{uJ,h) 
with that given in equation (|23p . The interaction term 
Ut{h) expressed in terms of the two other renormalised 
parameters, A(ft) and fj{h), is 



Utih) 



Km - 1) 

m{h) 



■Kh{fi{h) - 1) 



tan-i(^(^(M/A(ft)) 



(41) 



Detailed comparison of the RPT results for the transverse 
susceptibility based equations PO)) (or equivalently (US])) 
and (|^T|) with a direct NRG calculations have been given 
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earlier^^. They are in remarkably good agreement with 
the NRG results for all values of the magnetic field over a 
magnetic field range — fOTx < h < fOTx, asymptotically 
exact as a; — >■ 0, and satisfy the Korringa-Shiba relation. 

We can derive an approximation for the self-energy 
S:? (w, h) analogous to that given in equation (|35p . 



j:^iio,h)^U?{h) / Gf{iu + oj\h)xtiuj',h) 



duj' 



(42) 



Some RPT results based on equations ([^^ and (|40p for 
the self-energy and one-electron spectral density have 
also been given earlier I*' -'^^ and compared with the cor- 
responding results from direct NRG calculation. Again 
the agreement between the two sets of results was very 
good over a low frequency range — O.STk < uj < 0.5Tk 
for h = and over a larger range for higher values of h. 
These results demonstrate that it is possible to find an 
approximation using the RPT which will accurately re- 
produce the form of the self-energy in the low frequency 
regime for any value of the magnetic field. 



IV. RPT CALCULATIONS IN THE LOW FIELD 
LIMIT 

We need to extend the RPT calculations described in 
the previous section to see if, given with renormalised 
parameters for one field value, we can calculate the self- 
energy for another field value, and hence deduce how the 
renormalised parameters change as we vary the magnetic 
field. Before we consider this problem in detail we con- 
sider the simpler case of using the renormalised param- 
eters for ft, = to calculate the renormalised self-energy 
in the presence of a weak magnetic field h. As we want 
to keep the magnetic field term explicitly in the Hamilto- 
nian, and as the self-energy has the form ^^('^ + ah,h), 
we modify our procedure and use a slightly different form 
for the remainder term I]Jj°™(w + ah, h) via 



T.„{uj + ah.h) = I]^(0,0) 

d'E^{uj + ah,h) 
(UJ + ah) — 



j=h=0 



(43) 

S™"(w + aft,ft). 



As the renormalised parameters are those defined for 
h — 0, the renormalised self-energy and 4-vertex satis- 
fies the equations in pT|) and ([T2|) with h set equal to 
zero. Substituting this form into the equation ^ for the 
impurity Green's functions gives 



GAij,h) = 



LJ ~ id + ah + iA — SJ.(a;, h) 



(44) 



so the renormalised parameters now correspond to ft, = 0, 
and the renormalised self-energy is defined by 



"(a;,ft) = z|e^(w + CTft,ft) - i;<,(0,0) 



{uj + ah) 



dY,„{uj + aft, ft) 



Note that this self-energy for ft 7^ is different from 
the self-energy 'S{u!,h) defined earlier in equation ([5]), 
due to the different remainder term, so we use a slightly 
different notation to distinguish them. The renormalised 
self-energy S''(w, ft) is required to satisfy the equation in 
(ini) for ft = 0. 

We now show that exact result for E''(0,ft) to first 
order in ft is given by the self-consistent evaluation of the 
tadpole diagram as in mean field theory. The counter 
term Ai only cancels off the tadpole diagram for ft = 0, 
so taking this into account, we get for the renormalised 
self-energy, T,'^{uj, ft) = — ({/ + \z)m{h), with 



1 _j h~mQ,h) 

m(h) = —tan i 



(46) 



This equation for m{h) corresponds to the exact form 
given by the Friedel sum rule. Solving these mean field 
equations to first order in ft. we find for the zero field 
susceptibility. 



ttA 



XI 



2(1- (C/ + A3)AA) 



(47) 



This result for xi corresponds to the exact result given 
in equation (fT8| if f/ -|- A3 = Us- For ft = we also 
have Ut = Us so this mean field equation corresponds to 
the exact one given earlier in equation (j27p , but with the 
renormalised parameters A (ft) and Ut{h) taken at zero 
field. As the corrections to the zero field values for A(ft) 
and [/((ft) are of order ft^, it follows that equation ([77|) 
with the zero field renormalised parameters is exact to 
first order in ft. 

We can extend the calculations to include an uj- 
dependence in S^(a;, ft) by using the RPT equation 
with the zero field parameters in the transverse dy- 
namic susceptibility calculated from the NRG. From the 
change in the linear w-dependence of the renormalised 
self-energy with the magnetic field ft, we can calculate 
the change Az(ft) = z{h) — zifi). We can then compare 
the results with the corresponding value of Az(ft) de- 
duced from a direct NRG calculation of the self-energy. 
A comparison of the two sets of results is shown in fig- 
ure [5] for the case U/ttA — 3, ttA = 0.1. It can be seen 
that the two sets of results are in complete agreement in 
the very low field regime. This implies that, given the 
renormalised parameters at ft = 0, we should be able to 
calculate accurately both the change in z(h) and the shift 
in the real part of the self-energy with ft, provided we use 
a sufficiently small value of ft. This is all the information 
we require to calculate the renormalised parameters for 
small ft from the given values at ft = 0. 
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FIG. 6: A plot of the change in z{h), Az{h) — z{h) — z{0), as 
a function of /i/Tk calculated using the RPT for U/nA = 3, 
ttA = 0.1, with renormalised prameters at h — deduced 
from the NRG, compared with Az{h) deduced from a direct 
NRG calculation of the self-energy. 



V. SCALING EQUATIONS 

We now consider how to set up a scaling equation to 
deduce the renormalised parameters for a magnetic field 
h — Sh, given a renormalised self-energy that has been 
calculated using the renormalised parameters, fj{h), A(/i) 
and U{h) for a magnetic field value h. We first of all in- 
troduce a renormalised self-energy ^^(a;, h, Sh) following 
the steps from equation to (H5t but retain only the 
Sh term explicitly in the Green's function and the remain- 
ing part is absorbed to give the renormalised parameters 
fi{h) and A{h) as defined earlier, 

S''(a;, h, 6h) = z(/i){l]<,(cj, h-6h)- S^(0, h) - 



{uj — aSh)- 



du) 



=Sh-- 



J. (48) 



The self-energy S'^(aJ, h,Sh) satisfies the renormalisation 
conditions in ([TT]) for 5h — 0. The local Green's function 
for the system with a magnetic field h — 5h, using the 
parameters for the system in a field /i, takes the form. 



G„{uj,h^5h) = 



(49) 



z{h) 



w - aSh + ah{h) + iA{h) - S^(tJ, h, 5h) ' 



where h(h) = hf}{h). However, in terms of the renor- 
malised parameters for the field value h — Sh, this Green's 
function is 



Ga{u,h-6h) 



(50) 



z{h-Sh) 



UJ + ah{h - Sh) + iA{h - Sh) - Y.„{lu, h - Sh) 

We can find an expression for the renormalised parame- 
ters h(h — Sh) and z{h — Sh) in terms of the parameters 



for Sh = and the self-energy E^(w, /i, by equating 
the inverses of the Green's functions in equations (|49|) 
and (|50|) . Differentiating these with respect to oj, and 
putting w = 0, gives 



z{h~Sh) ^ z{h,Sh)z{h), 

where z{h,Sh) = 1/(1 - (0, /i, (5/i)). 
relation between A(/i — Sh) and A(/i), 



(51) 

Hence we get a 



K{h-Sh)^z{h,5h)K{h). (52) 

Equating the inverses of and (150]) and putting a; = 0, 
we find a relation between h{h — Sh) and h{K)^ 

h{h - Sh) = z{h, Sh){h{h) -Sh- E!f(0, h, Sh)). (53) 

We can find a alternative form of this equation by ex- 
panding Ti'^{0, h, Sh) to first order in Sh and then using 
the relation, 

d^{0,h,Sh)) 



dSh 



Sh=0 



^-p{0,h)U{h), 



(54) 



which follows from the Ward identity^i^SiiSi. We then get 
an equation for h{h — Sh) in the form, 

h{h - Sh) = z{h, Sh){h{h) - Sh{l + R{h))) + 0[{Sh)% 

(55) 

where R{h) is the Wilson ratio given by R{h) = 1-1- 
~p{0,h)tj{h). 

We also need to be able to calculate the new interaction 
term Ut{h — Sh). It can be seen from equation (PT|) that 
a knowledge of fi{h) — h{h)/h and A{h) is sufficient to 
determine Ut{h). Hence from the results for A(/i — Sh) 
and {h — Sh) we can deduce Ut{h — Sh). 

Finally, the full renormalised vertex U{h — Sh) can be 
calculated from the results for A(/i — Sh) and fj{h — Sh), 
by differentiating the magnetisation, as given in equation 
P^ . to determine the static longitudinal susceptibility, 
and then equating it to the expression given in equation 

m- 

We now have to amend the procedure outlined in the 
previous section to calculate the self-energy in the re- 
duced field S^(w, h, Sh) so that we can exploit equations 
([5^ and ([55]) or (|55p to extend the calculation of the 
renormalised parameters to smaller values of the field. 



VI. EXTENSION TO LOWER MAGNETIC 
FIELD VALUES 

The Green's function for a system with a magnetic 
field h — Sh, in terms of renormalised parameters for a 
field value h, takes a form similar to that given in equa- 
tion (j36p when we include the counter terms Xi^a{h) and 

A2,<t(Mi 



G„{iu,h-Sh) 



(56) 



z{h) 



uj + (T(h- Sh) + iA{h) + Ai,^ + ujX2^a - S;;(w, h, Sh) 



z(h)RPT 




-4 -2 2 4 

ln(h/T^) 



FIG. 7: A comparison the renormalised parameters, z{h) — 
A{h)/A, Ut(h)/U and fi{h) as calculated from the RPT and 
NRG for the symmetric Anderson model, with vrA = 0.1, 
U/tt/S. = 3 as a function of the logarithm of the magnetic 
field /i/Tk, where Tk = 7rA(0)/4. 
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FIG. 8: A comparison of the renormalised parameters, 
A{h)/A, U{h)/U and T){h) calculated using the RPT as a 
function of the logarithm of the magnetic field h/Tn, com- 
pared with the corresponding results from the NRG given in 
figure [1] 



The self-energy /i, ^/i) is now calculated with the 

quasiparticle propagator, 



Ga{u:,h^5h) 



(57) 



cj + (t(/i - 5h) + iA{h) + Xi,a + i^\2,a - h, 5h) 

The counter terms Ai.o- and A2,o-, however, are still de- 
termined by the conditions given in equations p7p and 
(|38| but with 5h = 0. 

We consider first of all the tadpole diagram, which 
now gives a finite contribution because it is not cancelled 
completely by the counter term for Sh 0. This di- 
agram, when the cancellation due to the counter term 
Ai^o- has been taken into account gives T,'^(uj,h,dh) — 
—aUt{h)ni{h,6h) where 



1 _^ -Sh-i:^o,h,dh) 

m[n,dn) — —tan i 

TT \ A 



(58) 



The mean field self-consistent equation for m{h, Sh) is 
then solved, and this term absorbed into the propagator 
in equation (j57p so the perturbation expansion is now 
about this mean field solution. The next step is to take 
the repeated quasiparticle scattering term into account 
in the calculation of the self-energy T:^{uj,h,5h) as in 
equation (|^^ . This constitutes our approximation for 
calculating Yi^{uj,h,dh). The renormalised parameters 
as a function of magnetic field can now be calculated for 
a field h — Sh using this approximation for Yi^{uj,h,Sh) 
together with the scaling equations ([5^ and ([55)) or (155)1 . 



0.5 
0.4 
0.3- 
0.2- 
0.1 


-0.1 



z(h) RPT 
U(h)/U RPT 
z(h) NRG 
C(h)/U NRG 



I 

// 



In(hn'j^) 







FIG. 9: The NRG and RPT resuhs for A{h)/A and f/(/i)/C/ 
as in figure [7] over the the low field regime. 



In figure[7|we compare the results for f]{h), Ut{h) and 
z{h) = A(ft,)/A obtained in the RPT, using the scal- 
ing equation to extend to the low field regime, with the 
NRG results. The results are for the symmetric model 
with [//ttA = 3, ttA = 0.1. It can be seen that the agree- 
ment with the NRG results is very good over the whole 
magnetic field range, especially the results for Ut{h) and 
z{h) = A(/i)/A. The results for fj(h) are smaller than 
the NRG results at lower field values but the difference 
is relatively small. In the very low field regime /i — >■ 0, 
the NRG results approach the value 2 corresponding to 
the Wilson ratio in the Kondo limit. The RPT results in 
this regime are smaller by approximately 3%. 

Within the approximation we have used we can cal- 
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FIG. 10: The renormalised parameters, ziK) = A(/i)/A, 
C/(/i)/7rA as calculated from the RPT for the same param- 
eter set as in figure [Jl showing to a good approximation the 
approach to the zero field strong coupling result U = ttA. 




FIG. 11: A comparison of I/ttAxs for /i = as deduced from 
the RPT (full curve) with the Bethe ansatz results (dashed 
curve) for a range of values of fZ/vrA. In the strong correlation 
regime I/ttAxs ATk/t^A. 



culate U{h), given Ut{h) and fj{h), rather more simply 
than the method described earlier, by summing the renor- 
malised equivalent of the diagrams in figure 2] to give 



Uih) 



Ut{h) 



i-Ut{h)nt{o,h,/^{h)) 



(59) 



Using the fact that Ilf (0, ft,, A(/i)) = mo{h,A{h))/h, we 
find U{h) = hUt{h)/h = ij{h)Ut{h). 

In figure [5] we compare the results for fi{h), U{h) and 
z{h) = A(/i)/A obtained in the RPT using the scaling 
with the NRG results as given in figure [T] for the same 
parameter set. There is excellent agreement between the 
U{h) from the RPT with the NRG result. As the values 



oizih) = A(/i) /A and U{h) are rather small in the strong 
coupling regime /i — >■ 0, we give an enlarged picture of 
the comparison of the results in the weak field regime 
in figure [3 The agreement of the RPT with the NRG 
results can be seen to be maintained down to values of 
the magnetic field several orders of magnitude less than 
the Kondo temperature Tk. Finally in figure [TUl we give 
a plot of the RPT values of U{h)/TTA and A{h)/A. The 
fact that the two curves merge as /i ^ shows that the 
RPT results asymptotically satisfy the relation U{0) = 
7rA(0), as given in equation (|28p corresponding to a single 
renormalised energy scale. 

In figure [TT] we compare the results for 1 /ttAxs from 
the RPT calculation in the zero field limit for a range 
of values of U/ttA with the corresponding Bethe ansatz 
results^^ given by 



1 



ttAxs (l + /)7rA' 



where 



and 



1/2 



^-ttU/BA+-!tA/2U 



I 



dx. 



(60) 



(61) 



(62) 



For U / vrA > 2, the integral term / is very small compared 
to unity so that in this regime I/ttAxs corresponds to 
4rK/7rA. 

It can be seen from this comparison there is very good 
agreement with the Bethe ansatz results in the strong 
correlation regime U /ttA > 2.5. It shows clearly that the 
RPT results give the correct form for Tk in the Kondo 
regime, not only in the exponential dependence on U but 
also in the prefactor. There is a progressive improvement 
in the agreement with the exact Bethe ansatz results with 
increase of U/ttA over the range 2 < U/ttA < 4.5. 

It might seem surprising to be able to get such pre- 
cise agreement by taking what would appear to be only 
a subclass of diagrams into account. The explanation is 
that we are considering the low energy regime and, for 
strong correlation, the dominant low energy scattering is 
with the spin fluctuations. The low energy charge fluctu- 
ations are suppressed in this regime and their effects can 
be taken into account by suitably renormalised vertices. 
For smaller values oil// ttA the charge fluctuations begin 
to play a role and have to be considered explicitly. In 
this regime they cannot be taken into account simply by 
the use of an effective frequency independent vertices. 



We are now in a position to calculate the low energy 
dynamic susceptibility and self-energy at ft, = and low 
fields entirely from the RPT. In figures [T^] and [T3] we 
compare the RPT results for the imaginary part of the 
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FIG. 12: A plot of lmxt(t^,0) versus lj/Tk for U/tvA = 3, 
ttA = 0.1 and h — 0, as calculated directly from the NRG 
(dashed curve) and from the RPT with renormalised param- 
eters also calculated using the RPT. 
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FIG. 13: A plot of Imxt(t^,/i) versus cj/Tk for J7/7rA = 3, 
ttA = 0.1 and h — Tk ~ 0.002, as calculated directly from 
the NRG (dashed curve) and from the RPT with renormalised 
parameters also calculated using the RPT. 



dynamic transverse susceptibility (full curve) with the 
corresponding results calculated directly from the NRG 
(dashed curve) for the case U/ttA = 3 and ttA = 0.1 and 
h = and h — Tk respectively. These have been cal- 
culated using equation (|23p for Xt('^j^) with the renor- 
malised parameters. It can be seen that they are in good 
agreement over the whole low frequency regime. The 
discrepancy in the peak height in figure [13] is almost cer- 
tainly due to the logarithmic broadening used in the NRG 
calculations which has the effect of reducing the height 
of any peak displaced from the origin. The higher the 
magnetic field value the larger the effect becomes as the 
peak gets shifted further from the origirJ^. 



In figures [TJ] and [TS] we make a similar comparison 
of the RPT results for the imaginary part of the self- 
energy calculated to second order in U for the same set 
of parameters, U/ttA — 3 and ttA = 0.1 and h ~ and 
h — Tk- The imaginary part of the self-energy T,^{uj, h) 
is related to the imaginary part of the renormalised self- 
energy Ecr(w,/i) from equation ([8|) via 

lmSa-{uj, h) ~ Imiig(aj, h). (63) 
z[n) 

There is very good agreement between the two sets of 
curves over the range |w| < 0.5Tk. The small discrepancy 
in the low frequency regime in figure [TS] could be due to 
the imaginary part of the self-energy in the NRG results 
does not precisely equal zero at w = 0. Surprisingly, using 
the RPT only to second order mU ioi h — extends the 
range of the agreement to < Tk. 

The low order RPT calculations give the asymptot- 
ically exact results corresponding to Fermi liquid the- 
ory in the low frequency regime. However, once the 
renormalised parameters have been calculated the renor- 
malised perturbation theory is completely defined, and 
can be extended to higher frequency and higher energy 
scales by taking higher order diagrams into account. 
The extension to the higher frequency range is currently 
being studied, and some preliminary results have been 
published^' 
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FIG. 14: A plot of the imaginary part of the self-energy 
E-f-(a;,0) versus cj/Tr for [//ttA = 3, ttA = 0.1 and /i = 0, 
as calculated directly from the NRG (dashed curve) and from 
the RPT with renormalised parameters also calculated using 
the RPT. 



VII. CONCLUSIONS 

We have set out to see whether we can access the low 
energy behaviour of a strongly correlated system using 
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FIG. 15: A plot of the imaginary part of the self-energy 
'Sf{Lo,h) versus lj/Tk for U/ttA = 3, ttA = 0.1 and h = Tk, 
as calculated directly from the NRG (dashed curve) and from 
the RPT with renormalised parameters also calculated using 
the RPT. 



a renormalised form of perturbation theory which would 
be applicable to a general class of models. We have taken 
as a test case the particle-hole symmetric single impurity 
Anderson model with parameters in the strongly corre- 
lated Kondo regime, where we have exact results which 
can be used to test any approximation used. The ap- 
proach has been based on the renormalised perturbation 
theory, which has previously been shown to give asymp- 
totically exact results for this model in the low energy 
regime in terms of renormalised parameters. Hitherto 
these renormalised parameters have been calculated in 
terms of the bare parameters that specify the model from 
an analysis of the low energy fixed point of an NRG 
calculation. The renormalised parameters can be deter- 
mined very accurately using this indirect approach, but 
this makes the method dependent on the NRG. Our goal 
has been to find an alternative and more general way to 
calculate these parameters, and to test the approxima- 
tions used by comparing the results with those derived 
using the NRG. To do this we have exploited the fact that 
the low energy spin fluctuations, which cause the strong 
renormalisation effects, can be suppressed by the appli- 
cation of a very strong magnetic field, allowing standard 
perturbation theory to be applied and the renormalised 
parameters to be calculated in the large field regime. The 
renormalised perturbation theory is then used to calcu- 
late the renormalised parameters on reducing the field 
/i by a small amount to h — 5h, so setting up a scaling 
relation. The solution of this scaling equation allows the 
field value to be extended down to ft, = 0. The approxi- 
mation used has been based on a mean field-like theory 
at each stage with a self-energy that includes the RPA- 
like fluctuations about the mean field. The results for 
the renormalised parameters have been compared with 
those calculated in previous work using the NRG for the 



complete magnetic field range. This relatively simple ap- 
proximation scheme gives remarkably good results for 
magnetic field values down to /i = 0, as demonstrated 
in the comparison with the Bethe ansatz results for the 
zero temperature susceptibility. This calculation demon- 
strates that it is possible to access the strong correlation 
regime using a perturbational analysis. 

There are some similarities with the approach used 
here for the Anderson model with two other perturba- 
tional approaches. In the local moment approach22i^ 
RPA diagrams are included in a two self-energy formal- 
ism, where the self-energies are constructed from a mean 
field broken symmetry state for ft. = 0. The symmetry, 
which should be restored by the low energy dynamics, is 
then imposed by the requirement that the imaginary part 
of self-energy is zero at w = 0, which determines the mean 
field magnetisation. This approach has an advantage in 
giving an interpolation such that the high energy fea- 
tures, corresponding to the atomic states, are included, 
and does give a low energy scale Tk which depends expo- 
nentially on U /ttA with an exponent in agreement with 
the Bethe ansatz result^. The fact that the symmetry 
has to be imposed, however, rather than develop natu- 
rally as the field is reduced in the RPT approach, is a 
limitation. 

The other related approaches have been based on the 
functional renormalisation group. A recent example 
is the work of Bartosch et alj^ who use a Hubbard- 
Stratonovich transformation and include both longitu- 
dinal and transverse spin fluctuations in calculating the 
self-energy. A cut off A was imposed to suppress the low 
energy spin fluctuations, which plays a role similar to the 
large applied fleld used in the RPT. Flow equations for 
the renormalised vertices as a function of A were then 
derived based using various approximations based on the 
truncation of the higher order vertices. The results, how- 
ever, for the quasiparticle weight factor z, which should 
behave as z ~ Tk/A in the Kondo regime, did not have 
the exponential dependence on [//ttA. 

The demonstration of the feasibility of this RPT ap- 
proach has been for a specific model, for a particular 
parameter set. The method, however, is a rather gen- 
eral one, and should be applicable to more general im- 
purity models. A good test case would be the general 
n-chaimel Anderson model with a Hund's rule exchange 
term. The renormalised parameters have been calculated 
for this model using the NRG for n = 32^, but for n > 2 
the NRG calculations become prohibitively difficult due 
to the sizes of the matrices to be diagonalised. It should 
also be possible to apply the approach to lattice models of 
strongly correlated systems, where in an appropriate ap- 
plied field the low energy behaviour of the system evolves 
continuously on reducing an appropriate applied field to 
zero. The next step is to test the approach for these more 
general types of models. 
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